Origin of spontaneous broken mirror symmetry of vortex lattices in Nb 
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Combining the microscopic Eilenberger theory with the first principles band calculation, we in- 
vestigate the stable flux line lattice (FLL) for a field applied to the four-fold axis; H \\ [001] in cubic 
Nb. The observed FLL transformation along Hc2 is almost perfectly explained without adjustable 
parameter, including the tilted square, scalene triangle with broken mirror symmetry, and isosceles 
triangle lattices upon increasing T. We construct a minimum Fermi surface model to understand 
those morphologies, in particular the stability of the scalene triangle lattice attributed to the lack 
of the mirror symmetry about the Fermi velocity maximum direction in k-space. 

PACS numbers: 74.25.Uv, 74.25.0p, 74.20. Pq 



A series of elemental metal superconductors, Nb, V, 
Al, • • • have a long history since its discovery of Hg in 
1911. Those are now regarded as a "conventional" su- 
perconductor where in fact the energy gap is isotropic 
and electron-phonon mechanism is known to work well. 
Thus it seems that no mystery remains there and every 
superconducting aspect is well understood and controlled 
since Nb is known to be wide practical applications from 
Nb-SQUID magnetometers to superconducting cavities 
for particle accelerators [ij. However, recent small angle 
neutron scattering (SANS) experiments on Nb discover 
plethora of the mysterious vortex morphology re- 
quiring us a consistent understanding of flux line lattice 
symmetry changes upon varying the applied field direc- 
tion continuously, that should satisfy a global constraint 
imposed by the hairy ball theorem[5|, namely any local 
change of FLL symmetry must be consistent with the 
change viewed globally. 

Among them recent finding by Laver et al. 0] is par- 
ticularly intriguing: They find a novel FLL whose half 
unit cell is a scalene triangle by SANS experiment on 
Nb under a field applied along H || [001] direction. By 
varying H and T they have succeeded in constructing 
the vortex lattice phase diagram (see Fig. 1(b)). Along 
the Hc2 line from low T and high H the square lattice 
tilted by 10.9° from the [100] direction changes into a 
scalene lattice followed by an isosceles lattice that is sta- 
bilized at higher T. Although the tilted square lattice has 
been known for a quite long time experimentally without 
any theoretical understanding Q , the discovery of the sca- 
lene lattice is new and intriguing since the scalene lattice 
breaks the fundamental mirror symmetry spontaneously 
in cubic crystalline Nb. Notice that among various su- 
perconductors, conventional and unconventional, there is 
no known case to exhibit either scalene triangle or tilted 
square FLL under the four-fold symmetric field config- 
uration, namely, tetragonal crystals under H [| [001] 
such as CeCoIns 0, Sr2Ru04 [1], La1.83Sro.12O4 0, 
ErNiaBaC 0, YNizBaC [HI, and TmNiaBaC [ll]. 

According to Takanaka [13] who is a pioneer in this 
field of theoretical FFL morphology studies, the free en- 



ergy of FLL contains the so-called phasing energy that 
depends on the relative orientation between the FLL and 
underlying crystal lattice. This implies that the observed 
tilted square lattice may be stabilized by this phasing en- 
ergy. This is expanded in terms of the higher order har- 
monics of the four-fold symmetry {9 is the angle from the 
[100] direction), Fphasmg = S„=i,2, --44n cos4n6' whose 
minima may give rise to the tilting angle 9 = 10.9° from 
the [100] direction, meaning that at least two or more 
higher harmonics are simultaneously no n- vanishing. 

Here the purposes of this paper are two-fold: The first 
principles band calculation based on density functional 
theory (DFT) yields precise information on the electronic 
properties of the normal state in general. By combining 
it with the microscopic Eilenberger theory, we can estab- 
lish a truly first principles framework for a superconduc- 
tor, in particularly for vortex state under an applied field 
without any adjustable parameters. In order to bench- 
mark it we apply this to explain those intriguing vortex 
morphologies, which will turn out that this framework is 
remarkably successful. 

Thus the second purpose is to understand the phys- 
ical origin for the scalene FLL structure by construct- 
ing a minimum model Hamiltonian to describe the FLL 
transformation mentioned above. The mirror symmetry 
breaking with the scalene FLL is attributed to the spe- 
cial character of the Nb electronic band structure. This 
study points to the fact that our combined framework 
of DFT and the Eilenberger theory paves the way to 
truly first principles theory for superconductors proper- 
ties under a field without adjustable parameter. Previ- 
ously several attempts have been done to perform this 
pro gram , but those are limited to either discussions on 
iJ^9[14j] or resorting to approximate solutions for Eilen- 
berger equation [lit. 

In order that our calculations are tractable within a 
reasonable computational time frame, we limit our dis- 
cussion along Hc2, that reduces significantly the compu- 
tational time because we must minimize the free energy 
in the multi-dimensional parameter space spanned by the 
FLL unit cell shape and its orientation relative to the un- 
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der lying crystal axes. 

The quasiclassical Eilenberger theory is quanti- 

tatively valid when ^ ^ l/fcp (^f is the Fermi wave num- 
ber and ^ is the superconducting coherence length). The 
quasiclassical Green's functions (7(w„,k, r), /(cj„,k, r), 
and /'''(wn, k, r) are calculated by the Eilenberger equa- 
tion 

{w„ +VF • (V + iA)}/(w„,k,r) = A.g(cj„, k, r), 
{w„ ^VF-{V- iA)} /tK, k, r) = A*5K, k, r),(l) 

where g — (1 — //^)^^^, Reg > 0, A is the vector potential 
and Vi? is the Fermi velocity. At t ~ T/Tc, the self- 
consistent gap equation is given by 



Alnt = i^(^(/K,k,r))k-^) 



(2) 



where (• • •)k = / deD{e, ej.)(- • •) with DiO, ep) being the 
angle-resolved density of states (AR-DOS) at the Fermi 
surface. 

Along Hc2, we linearize the Eilenberger equation as 



VF-(V + iA)]/K,k,r) = A. 



(3) 



By expanding A and f{ujm k, r) in terms of the Landau- 
Bloch function V'jv.qlr) whose 7V-th coefficients are Ajv 
and /^{ujn, k), we can obtain the eigenvalue equation for 
A, from which Hc2, Ajv and /A'(a;„, k) are evaluated [20l|. 
Those eigenfunctions lead to the gap function A(r) and 
the quasiclassical Green's function /(a;„,k, r). 

In order to compare various FLL forms to find the most 
stable vortex configuration near i?c2, we have to calculate 
the free energy F given by 



F = k^V X A)2 + |A|2 1ni- 



|A|2 



-{I) 



(4) 



with I = {g - l)[2w„ + vf ■ (Vln(///t) + 2iA)] + 
(/A* -I- A/t). Via the Abrikosov identity and expand- 
ing the normalization condition up to the next order; 
g ^ 1 — — (//^)^/8, we obtain the free energy valid 

near Hc2 as Il8| 



F 



F = 



i?2 
t 



4^2 /l2 ^ 



F + 1 



(5) 



hi 



A/t))k-^, (6) 



where (• • •) denotes the spatial average within a unit cell, 
and hs is the magnetic field induced by supercurrent. 
B — H+hs and k'^ = 7C(3)k2/18. Since we know that the 
vortex morphology is independent of GL parameter k 21| , 
here we have used a large k. The free energy minimum 
corresponds to the minimum of F that gives the stable 
vortex lattice configuration. 





FIG. 1: (Color online) Calculated Fermi surfaces (a) and 
(b) for Nb in first Brillouin zone where the former (latter) is 
closed (extended). 




FIG. 2: (Color online) The cross sections perpendicular to 
the field direction H|| [001] for the Fermi velocity anisotropy 
vf{0) (a), AR-DOS D{e,eF) (b) and the inverse of the AR- 
DOS (c) where the red line drawn on its maximum direction 
indicates the lack of the mirror symmetry about this a^ds 
tilted by 11° from [100]. 



The Fermi velocity -VFiO) and AR-DOS D{e,eF) are 
evaluated by using tetrahedron method within DFT^^. 
The first Brillouin zone is divided into 135'^ parallel pipes, 
each of which is further subdivided into 6 pieces. As 
shown in Fig. [U our DFT calculation reproduces the 

Then the band in- 



well-known two Fermi surfaces 22 



formation is used to yield the Fermi velocity anisotropy 
vf{0) shown in Fig. ^a.) and the AR-DOS D{e,eF) in 
Fig. m^b) where we map the three dimensional objects 
onto the cross sections perpendicular to the field direc- 
tion for the study of vortex states. Those quantities are 
used for evaluating the free energy. Here we also show 
the inverse of AR-DOS in Fig. ^c) for later purpose. 

Figures [3] and [4] show the FLL transformation along 
Hc2 where the unit cell shape characterized by the an- 
gles, a, /3 and 7 and the tilting angle 5 from [100] are 
displayed in Fig. 3(a) and compared with the observed 
phase diagram^] in 77 vs T for Nb in Fig. 3 (b). Note 
that the square lattice tilted by 5=14° is stabilized at the 
low T and high H region. This tilting angle nearly coin- 
cides with the observed value 10.9° and their H—T region 
also coincides with the observation. Upon increasing T, 
the tilted square lattice changes into the scalene lattice 
at < = 0.3. Then after gradual transformation, keep- 
ing the scalene nature intact for a finite H region. This 
scalene lattice that breaks the mirror symmetry of FLL 
finally gives way to the isosceles triangle lattice at higher 
T. This isosceles is orientated along the [100] direction 
(5=0°) that coincides with the observation. This lock- in 
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FIG. 3: (Color online) Calculated FLL transformation along 
Hc2 (a) where the square lattice is tilted by 5=14° from [100], 
the scalene FLL is stabilized in intermediate T, and the isosce- 
les triangle in high T is along [100]. Observed phase diagram 
in Nb (b). 
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FIG. 4: (Color online) Detailed FLL transformation of Fig.l 
(a) low T from tilted square to scalene triangle and (b) inter- 
mediate T from the scalene to isosceles triangle. The insets 
(c) and (d) show the experimental dataf^] for reference: (c) 
data at 2K as a function of H, and (d) data at _H'=200mT as 
a function of T. 



transition occurs smoothly. Upon further increasing T, 
this isosceles triangle lattice tends to the equilateral tri- 
angle toward Tc smoothly. The overall FFL transforma- 
tion characterized by the unit cell shape and orientation 
relative to the crystal lattice along Hc2 almost perfectly 
reproduces the observation0,|3]. It should be emphasized 
that there is no adjustable parameter here and all infor- 
mation on the electronic structures in the normal state, 
including the Fermi velocity anisotropy and the AR-DOS 
is coming from the first principles DFT calculations. 

In order to understand the physical origins of the re- 
markable success of the above Eilenberger theory com- 
bined with DFT, we construct a minimum model for 
grasping the essential features of those results. Namely 
(1) the tilted square at low T, (2) the scalene triangle in 
the intermediate T and (3) the isosceles triangle whose 
nearest neighbor is along [100] at higher T. We model 
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FIG. 5: (Color online) FLL transformation based on the 
minimum model with Pi — 0.1 and /3g — 0.3. The low T 
square lattice tilted with 18° from [100] axis transforms into 
the scalene FLL that is followed by the isosceles FLL at higher 
T where the nearest neighbor orients along [100] direction. 



the Fermi velocity anisotropy in the Fermi circle whose 
radius is kp, that is, vp — l3{6)vpU. (3{9) parametrizes 
the Fermi velocity anisotropy, and u is the unit vector 
on the plane perpendicular to the field direction. In this 
model, the AR-DOS is given by D{9,eF) = l/P{e)v%. 
Namely, the Fermi velocity anisotropy and AR-DOS are 
not independent quantities as in DFT. We can expand as 
I3{9) = S„=o,i,2, - /34n cos(4n6'). Since it will be clear that 
/34 only model cannot account the above features (l)-(3), 
the minimum model consists of jS^ and /^g- Namely, we 
consider the following model: 



vi.(0) =4u/?(0) 



i;^u(l-f;34Cos46i-/38( 



We have performed extensive calculations based on this 
model, plugged it into the Eilenberger equation above to 
find the conditions to reproduce the above features. Then 
we have found the necessary condition for scalene triangle 
to appear as /Jg > /54 > that also explains the features 
(1-3) simultaneously. 

Figure 5 shows an example for this case {Pa = 0.1 and 
/Sg = 0.3). In the lowest T region, the square lattice 
appears with the tilting angle = 18° while the Fermi 
velocity maximum is at = 21.3°. There is a rule that 
at lower T the stable square lattice is oriented along ei- 
ther the AR-DOS maximum for isot rop ic gap case or the 
nodal direction for d-wave gap case[2lj. 

Then by increasing 7J, the tilted square deforms con- 
tinuously into the scalene triangle FLL as seen from Fig. 
5 in the middle T. This scalene lattice is locked-in the 
isosceles triangle lattice oriented along [100]. Those se- 
quences of the FLL transformation are same as in DFT 
case above. It is now clear that the essential ingredi- 
ent for the scalene triangle to appear is the lack of the 
mirror symmetry about the Fermi velocity maximum di- 
rection as indicated in the inset of Fig. 5 and also note 
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(a) 




FIG. 6: (Color online) FLL transformations as T increases 
for the two cases where the bold curves are the Fermi velocity 
anisotropics: (a) FLL changes from the square to the scalene 
for the minimum model with /?4 = 0.1 and j3s, = 0.3. (b) 
FLL changes from the square to the isosceles for /34 = 0.1 
where the mirror symmetry is unbroken about the maximum 
direction [100]. 



that in Fig. 2(c) the inverse of AR-DOS also shows this. 
This velocity maximum simultaneously explains why the 
square lattice is tilted. Namely the tilted square FLL 
and scalene FLL are intimately tied with each others. 
Note that while our P^-jSg, minimum model explains the 
appearance of the scalene FLL, the stable T region is far 
off the observation. 

In order to understand the origin of the scalene FLL, 
we display the deformation of this lattice in Fig. 6(a). It 
is seen from it that the deformation occurs so that one 
of the two diagonal axes of the parallel piped unit cell 
continuously rotates, keeping the other axis intact. This 
movement can be interpreted as the scalene FLL seeks 
stable configuration by adjusting the unit cell angle and 
its orientation by utilizing the possible freedoms allowed. 
As a result of the broken mirror symmetry about the ve- 
locity maximum direction the scalene FLL could appear. 
This situation is quite different from that in the /34 only 
(the same for the /?§ only model) where the mirror sym- 
metry is preserved in the velocity anisotropy along the 
maximum direction. As shown in Fig. 6(b), the FLL de- 
formation from the square to isosceles triangle proceeds 
as follows: The two diagonal axes of the parallel piped 
unit cell keep unchanged, always orthogonal with each 
other and the two angles of the isosceles continuously 
changes. 

It is important to notice that one of the two diagonal 
axes is kept to point along the direction of the veloc- 



ity minimum during the deformation. In other words, 
the base of the isosceles is always kept along the veloc- 
ity minimum direction. The pairs of the vortices on the 
unit cell move along those symmetry constraint axes, pre- 
serving the mirror symmetry. This is contrasted with the 
present scalene case where the two diagonal axes rotate 
each other, thus no symmetry constraint during the de- 
formation because the mirror symmetry about the maxi- 
mum direction is broken from the outset. Comparing the 
two cases with either broken mirror symmetry or with un- 
broken symmetry, it is concluded that a necessary condi- 
tion, but not sufficient condition for the scalene triangle 
FLL to occur is that in the anisotropic velocity the mirror 
symmetry is broken about the anisotropy maximum di- 
rection. The broken mirror symmetry of the FLL that is 
real space is originated from the mirror symmetry break- 
ing in the reciprocal space. Furthermore, it is interesting 
to notice the experimental fact that the appearance of 
the scalene FLL is tied with the appearance of the tilted 
square lattice in Nb. As mentioned above, the square 
FLL is oriented along the velocity maximum direction 
at lower T, that is a general theoretical fact within our 
model calculations. Thus, in Nb the Fermi velocity max- 
imum is situated along the direction tilted by 11° from 
the [001] axis and the mirror symmetry about this axis 
must be broken. 

In summary, we demonstrate that the first principles 
band calculation combined with microscopic Eilenberger 
analysis yields an excellent explanation of the observed 
FLL transformation in Nb under H \\ [001], including 
the tilted square and scalene lattices without adjustable 
parameter. This opens a door that this combined method 
provides a fundamental framework to truly understand 
the vortex matter from first principles. We also show the 
possible reasons why the spontaneous mirror symmetry is 
broken in the scalene lattice by constructing a minimum 
model for the Fermi velocity anisotropy. 
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stay in Aspen Center for Physics. 
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